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Résumé 

Nous exposons systématiquement les relations entre les modèles 
de gaz avec particules dures, l'énumération des animaux dirigés, et 
l'algèbre des séries formelles en variables partiellement commutatives, 
d'après les travaux de X.G.Viennot. Nous donnons une solution com- 
plète simplifiée, en utilisant cette algèbre, de l'énumération des ani- 
maux dirigés sur réseau plan. L'article inclut un programme rapide de 
génération aléatoire de ces animaux, avec distribution uniforme, et des 
images obtenues à l'aide de ce programme. 

Abstract 

We give a systematic présentation of relations between lattice gas 
models with hard-core interactions, enumeration of directed-site ani- 
mais, and the algebra of formai power-series in the partially commuta- 
tive case, along the work of X.G.Viennot. We présent a complète and 
simplified solution, using this algebra, of the directed animais enumer- 
ation problem in two dimensions, including a fast program for random 
génération of such animais (with uniform distribution) and images pro- 
duced by this program. 
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1 Introduction 



Certains physiciens (voir par exemple 1^, ^, |9|, |ï^, [ï^) s'intéressent à 
des modèles thermodynamiques extrêmement schématiques, mais où. appa- 
raissent déjà des phénomènes critiques, dont ils peuvent faire une étude 
détaillée. Dans ces modèles, des particules sont placées sur certains sommets 
d'un graphe G, appelés sites (au plus une particule par site) ; la disposition 
des particules forme une configuration, et V énergie d'une configuration est la 
somme de deux séries de termes : les premiers représentent la contribution 
individuelle de chaque particule, les seconds V interaction entre particules 
voisines (c'est à dire reliées par une arête du graphe). 

Les configurations étudiées sont de très grande taille, et on s'intéresse à 
leurs propriétés statistiques ; cette étude se ramène à celle de la fonction de 
partition du modèle, notée traditionnellement Z, qui joue un rôle voisin d'une 
fonction génératrice en combinatoire ; d'ailleurs le calcul de Z se ramène 
souvent explicitement à celui d'une série énumératrice : voir en particulier 

m- 



X.G.Viennot a découvert - voir par exemple [|Î6[ - que, dans le cas 
particulier des modèles dits avec particules dures, cette étude est liée à 
celle du monoïde partiellement commutatif, introduit en 1969 par P.Cartier 
& D.Foata @, et qui fait l'objet aujourd'hui de recherches actives dans 
un domaine à priori éloigné, celui de la modélisation des algorithmes par- 
allèles en informatique théorique. Cet article expose cette approche de façon 
systématique, et inclut la solution selon cette ligne, due aux auteurs, du 
problème de l'énumération des animaux dirigés plans. 

La section ^ présente la classe de modèles étudiés, et contient les défini- 
tions issues de la physique, en particulier celle de la fonction de partition. 

La section ^ introduit la notion empilement de configurations, les opéra- 
tions algébriques sur ces objets, les théorèmes fondamentaux de cette algèbre 
(qui est celle des séries en variables partiellement commutatives), et leurs 
corollaires pour les modèles avec particules dures. 

La section |^ présente les animaux dirigés sur un réseau, et montre que, 
sous certaines conditions, l'énumération de ces animaux, sur un réseau de 
dimension d, se ramène directement à celle d'empilements, sur un réseau de 
dimension d — 1. 

Dans la section ^, on achève les calculs dans le cas le plus simple, celui du 
modèle de gaz sur réseau linéaire, qui correspond aux animaux dirigés sur 
réseau plan. Dans ce cas, les séries en variables partiellement commutatives 
vérifient des équations algébriques simples, qui fournissent à la fois les séries 
énumératrices des animaux, et des bijections avec des objets combinatoires 
classiques, les chemins (ou mots) de Motzkin. 

Ces bijections sont utilisées section ^, pour programmer la construc- 
tion d'animaux dirigés plans aléatoires, avec distribution uniforme, en temps 
linéaire. La section |6| est illustrée de dessins typiques d'animaux de grande 
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taille, produits par notre logiciel ; leur aspect filiforme surprend l'intuition. 

La section |^ énumère les principaux problèmes ouverts, en particulier 
celui d'une interprétation combinatoire de la solution de R. Baxter pour le 
modèle des hexagones durs. 

2 Modèles en mécanique statistique 

2.1 Modèle d'Ising et modèle de gaz 

Dans le modèle d'Ising du ferromagnétisme, chaque site i est occupé par 
un petit aimant, avec deux orientations opposées possibles, notés ai = ±1, 
et l'énergie est donnée par : 

E = -HY,a,-J Y. (1) 

site i arête 

OÙ. H désigne le champ magnétique externe, et J la constante d'interaction. 
L'énergie due au champ externe est minimale lorsque les aimants sont ori- 
entés dans le sens du champ ; lorsque la constante J est positive, l'énergie 
d'interaction est minimale lorsque tous les aimants sont orientés dans le 
même sens, d'oij le nom de ferromagnétisme — on ne s'étonnera pas de con- 
stater que l'énergie minimale est négative, car E est défini à une constante 
additive près. Lorsque J < 0, chaque aimant pousse ses voisins à s'orienter 
en sens contraire de lui-même : on parle d'antiferromagnétisme. 

Il existe un modèle analogue des gaz : chaque site i est occupé ou vide, 
ce qu'on note Si = 1 ou 0, et l'énergie vaut : 

-E = -/i ^ Si - e ^ SiSj (2) 

site i arête 

fi est appelé potentiel chimique ; le cas e > correspond à une attraction 
entre particules, et e < à une répulsion ; le cas e = est celui, sans intérêt 
ici, du gaz parfait. 

Dans le modèle du ferromagnétisme, il existe une interaction entre chaque 
paire de sites voisins ; dans le modèle du gaz, cette interaction n'existe que 
si les sites sont occupés. Pour comparer les deux modèles, on effectue le 
changement de variable ai = 2si — 1 ; l'énergie d'interaction associée à l'arête 
devient alors : 

fTjCTj = 4:SiSj — 2si — 2s j + 1 

Lorsqu'on somme cette quantité pour toutes les arêtes, le terme linéaire —2si 
apparaît avec un coefficient égal au degré qi du site i (le degré d'un sommet 
est défini comme le nombre d'arêtes incidentes). Donc (|ï]) devient : 

E = -2H Y Si + 2JY^ QiSi - 4J J2 ^i^j 

site i site i arête 
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en ignorant la constante nH — mJ, où n (resp. m) désigne le nombre de 
sommets (resp. arêtes) du graphe. Si le graphe est régulier, chaque degré qi 
est égal à une constante q, et on voit que dans ce cas les deux modèles sont 
équivalents, avec : 

I2 = 2H - 2qJ, e = 4J 

(cf. |, section 1.9]). 



2.2 Fonction de partition 

L'axiome essentiel de la mécanique statistique (formule de Gibbs), af- 
firme que la probabilité d'observation d'une configuration donnée, d'énergie 
totale E, est proportionnelle à e~^^, où (3 est une constante positive. /3 
dépend du milieu extérieur, dans lequel est plongé le modèle ; plus préci- 
sément : 

oii k désigne la constante de Boltzmann, et T la température absolue. Lorsque 
[3 est grand, les configurations de forte énergie deviennent très improbables : 
le milieu extérieur est froid ; inversement, lorsque /? tend vers zéro, toutes 
les configurations tendent à devenir équiprobables : le milieu extérieur est 
chaud. 

La fonction 

Z(T) = ^e-^^ = ^e-^/^^ 

où la sommation est étendue à toutes les configurations possibles, est appelée 
fonction de partition — définir l'énergie E à une constante additive près, a 
pour effet de définir Z à une constante multiplicative près. La probabilité 
d'observer une configuration donnée, d'énergie E, vaut donc : 



Z{T) 



(3) 



La température n'est pas en général la seule variable dont dépende la 
fonction Z : d'autres grandeurs physiques, comme le champ magnétique, 
peuvent intervenir dans le calcul de l'énergie d'une configuration, et donc 
de Z. Les physiciens montrent que l'étude statistique du modèle se ramène 
toujours à l'étude de la fonction Z et de ses dérivées ou dérivées partielles ; 
la section suivante illustrera ce principe général. 



2.3 Modèles avec particules dures 

Un cas limite des modèles de gaz précédents est obtenu en supposant que 
la répulsion entre particules est si forte que les configurations comportant 
des particules voisines ont une probabilité négligeable d'exister ; autrement 
dit, chaque particule est suffisamment dure, et occupe un volume tel, que 
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deux sites voisins ne peuvent être occupés simultanément. Les seules con- 
figurations possibles correspondent alors aux ensembles de sommets qu'on 
appelle stables en théorie des graphes : deux sommets distincts d'une même 
configuration ne sont jamais reliés par une arête. L'énergie d'une telle con- 
figuration de n éléments se réduit à — n/i, et en posant 



la fonction de partition vaut : 

Z(t) = ^a„r (4) 

oii a„ désigne le nombre de stables de taille n. En utilisant @, on en déduit 
que le nombre moyen de particules vaut : 

Unant^ tZ\t) d 

Cet exemple est caractéristique du rôle de la fonction de partition Z et 
de ses dérivées dans un calcul de moyenne. 



2.4 Limite thermodynamique 

La définition de la fonction de partition pour un graphe infini n'est pas 
immédiate : les formules précédentes n'ont plus de sens, puisque le nombre de 
configurations est infini. D'autre part la seule classe intéressante de graphes 
infinis, dans ce contexte, est celle des réseaux. 

En mathématiques (voir par exemple Q), un réseau de dimension d 
désigne habituellement, dans un espace affine, un ensemble de points iso- 
morphe à Z'^ ; les points associés, par l'isomorphisme, à ceux dont les coor- 
données valent ou 1, forment le parallélotope de base B ; dans un espace 
euclidien, un réseau est cristallographique si les symétries orthogonales par 
rapport aux faces de B laissent le réseau invariant. 

Dans le présent contexte, un réseau désigne un graphe, dont l'ensem- 
ble des sommets est un réseau cristallographique ; le choix le plus simple 
pour les arêtes est de ne relier que les voisins les plus proches (réseau nn, 
pour "nearest neighbour", chez les anglo-saxons), mais on peut ajouter des 
arêtes qui relient les voisins "au second rang" (réseau nnn, pour "next near- 
est neighbour"), etc. . . A un même réseau cristallographique peuvent donc 
correspondre plusieurs graphes : sur la figure ||-(a), chaque sommet est de 
degré 4, tandis que sur la figure 0-(b), chaque sommet est de degré 8. 

Notons qu'on peut étendre la notion de réseau en affaiblissant les hy- 
pothèses sur l'ensemble des sites : par exemple, si on pave le plan par des 
hexagones, les sommets ne forment plus un ensemble isomorphe à (il 
manque les centres des hexagones), mais le graphe obtenu est régulier, et 
définit un modèle d'interaction intéressant, oii chaque site est de degré 3. 
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(a) (b) 
FiG. 1: Réseaux carrés (a) nn (b) nnn 



Pour définir la fonction de partition d'un réseau G, on considère l'inter- 
section Gk de G avec un disque de rayon donné, oii k désigne le nombre de 
sites de Gk- Soit la fonction de partition de Gk ; Zj^ n'est pas une quantité 
proportionnelle à la taille k du système (le nombre de configurations croît 
exponentiellement avec la taille), mais la mécanique statistique enseigne que 
c'est le cas de InZk, lorsque k est grand; plus précisément, la fonction Z 
définie par : 

In Z = lim — — ^ (6) 

fc^oo k 

existe, et cette limite est indépendante de la définition précise de Gk- Les 
physiciens appellent ce type de limite la limite thermodynamique ; consulter 
R.Baxter [Q, ou un cours de mécanique statistique, comme Q, pour une 
présentation complète. 

Les phénomènes critiques correspondent aux singularités de la fonction 
de partition Z définie par (^), et à celles de ses dérivées ou dérivées partielles. 
Le nombre moyen de particules par site, appelé densité moyenne par site, 
est donné, sur un réseau infini, par le membre droit de la formule (|5|) oii Z 
désigne cette fois la fonction de partition pour la limite thermodynamique. 

Le modèle des carrés (resp. hexagones) durs correspond à l'étude d'un 
gaz, dont les particules sont disposées sur un réseau carré (resp. hexagonal) 
infini, de telle sorte que deux sites voisins ne peuvent être simultanément 
occupés. Pour rendre visible cette dernière condition, on peut représenter 
les particules par des carrés (resp. hexagones) qui ne se chevauchent pas, 
comme sur la figure ||. Dans la suite de cet article, un site occupé par une 
particule sera simplement indiqué par un disque noir. 

Le calcul de la fonction de partition du modèle des hexagones durs a été 
réussi par Baxter en 1980 : voir |^ pour une présentation complète ; c'est 
un calcul complexe, qui fait intervenir entre autres les identités de Rogers- 
Ramanujan ; le calcul de la fonction de partition du modèle des carrés durs 
reste un problème ouvert. 
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3 Séries en variables partiellement commutât ives 



3.1 Empilements 

Soit un graphe G ; les sommets de G formeront dans la section suivante 
les lettres d'un alphabet, et leur ensemble sera donc noté A ; une arête de G 
est une paire (non orientée) {a, b} de sommets distincts. 

Deux sommets de G sont dits voisins s'ils sont identiques, ou reliés par 
une arête de G ; le voisinage V{a) d'un sommet a est l'ensemble des som- 
mets voisins de a ; le voisinage V{B) d'une partie i? Ç ^ est défini comme 
d'habitude par : 

ViB) = U V{b) 

beB 

Dans toute la suite nous conviendrons qu'une configuration C est un 
ensemble de sommets tel qu'aucune paire de sommets distincts de G ne 
soient voisins (c'est ce qu'on appelle un stable en théorie des graphes). Un 
empilement est une suite finie (Ci, C2, . . . , C^) de configurations non vides 
telle que : 

Ci Ç F(a_i), l<i<n (7) 

L'entier n est la hauteur de l'empilement ; chaque configuration Gi est ap- 
pelée couche numéro i, et la première couche, Ci, est la base de l'empilement. 
Si n = 0, l'empilement est vide, et ne possède pas de base. 

Un empilement peut aussi être décrit comme un ensemble de cellules 
(a, i), où a est un sommet de C, et i un entier inférieur ou égal à la hauteur 
n de l'empilement ; i est la hauteur de la cellule, et la couche numéro i est 
formée des cellules de hauteur i ; les cellules ayant même projection o sur A 
forment la fibre Fa au-dessus de o. On dira pour abréger que deux cellules 
(a, i) et (6, j) sont voisines si les sommets a et 6 sont voisins; on peut alors 
reformuler la définition d'un empilement de la manière suivante : 

- deux cellules distinctes de la même couche ne sont jamais voisines (8) 

- pour i > 1, chaque cellule de la couche i possède au moins une voisine 
dans la couche i — 1 (9) 

La taille \E\ d'un empilement E est égale, par définition, au nombre de 
cellules qui composent E. Une configuration peut être confondue avec un 
empilement de hauteur 1 ; à la configuration vide, correspond l'empilement 
vide, de taille et de hauteur nulles. 

Exemple 1: la figure |3| montre un empilement de hauteur 4, et 
de taille 8 : 



Cl 


= {a,c} 


U(Ci) 


= {a, 


b,c,d,e,g] 


C2 


= {b,e,g] 


V{G2) 


= {a, 


b,c,e,f,g,h} 


C3 


= {e} 


V{C^) 


= {a, 


e,f,h} 


C4 


= {a,/} 
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g 






e 


f 









G 




Co 



Cl 



E 



FiG. 3: Un empilement E de configurations sur le graphe G 

Cl ) (^2 ) C'a ) C4 sont des configurations qui vérifient la condition 
(0), donc 

E = (Cl, C2, C3, C4) 
est un empilement, constitué des cellules : 

(a, 1), (c, 1), {b, 2), (e, 2), {g, 2), (e, 3), (a, 4), (/, 4). 



La notion d'empilement est due à X.G.Viennot ||Î5|, qui représente les 
cellules par des pièces géométriques, dessinées de telle sorte que deux pièces 
aient une intersection non vide si et seulement si elles sont voisines. Les 
pièces de la couche Ci sont posées sur le sol, celles de la couche C2 sont 
posées par dessus, et ainsi de suite (deux pièces de la même fibre ont même 
projection sur le sol). Grâce à la condition (^), chaque pièce d'un empilement 
repose sur sa (ou ses) voisine(s) de la couche inférieure. 



3.2 Monoïde partiellement commutatif 

Soit A* le monoïde libre sur l'alphabet A des sommets du graphe G ; le 
monoïde partiellement commutatif associé à G, noté pc(G), est le quotient 
de A* par les relations de commutation : 



ab = ba 



pour tous sommets a, b non voisins dans C (on dit que G est le graphe de 
non-commutation). Les éléments de pc(G) sont appelés traces, et la trace 
du mot u est notée [u]. Deux mots u,v de A* ont même trace si l'on peut 



11 



passer de l'un à l'autre par une suite de commutations autorisées, c'est à 
dire d'échanges entre lettres consécutives désignant des sommets de G non 
voisins ; dans ce cas on dit qu'ils sont équivalents, et on note : u ^ v. 

A une configuration C = {a, b,c, . . .}, associons la trace du mot u{C) = 
abc ... ; par définition d'une configuration, les lettres a,b,c, . . . commutent 
toutes entre elles, et la trace [ti(C)] est donc bien définie, car elle ne dépend 
pas de l'ordre dans lequel on énumère les éléments de C. 

A un empilement E = {Ci, C2, ■ ■ ■ , Cn), associons le produit des [n(Ci)] 
dans pc(G) (si l'empilement est vide, on lui associe la trace du mot vide), 
et notons <î> cette application. 

Exemple 1 (suite): la trace $(£') de l'empilement E de la fig- 
ure ^, est la classe du mot acbegeaf modulo les 16 commuta- 
tions : ac = ca, af = fa, . . . , fh = hf. 

Réciproquement, à tout mot u = ai ... an de j4*, associons l'empilement 
E{u) constitué des cellules (oj, h{ai)), oii la hauteur h{ai) de la lettre est 
définie par récurrence comme suit : 

h{ai) = 1 + max{ /i(aj) | j < i et aj G V{ai)} (10) 

avec la convention max(0) = 0. 

On vérifie facilement que E{u) est bien un empilement - c'est à dire 
satisfait aux conditions ^ et (^) ~ et que cet empilement ne change pas si 
l'on applique à u une commutation autorisée. Par passage au quotient, on 
obtient donc une application ^ de pc(G) dans l'ensemble des empilements, 
et on a le résultat suivant, énoncé dans ||ï^ : 

Théorème 1 Les applications ^et^ sont des bijections, inverses l'une de 
l'autre, entre l'ensemble des empilements sur un graphe G et le monoïde 
partiellement commutatif pc{G) . 

Preuve: lorsqu'on associe un mot n à un empilement E, en énumérant les 
couches de E, puis qu'on calcule les hauteurs des lettres de u par la formule 
(10), on retrouve les numéros des couches. Inversement, lorsqu'on part d'un 



mot u, la définition (|l^) implique que les lettres de hauteur 1 commutent 
avec toutes celles qui les précèdent, et donc u ^ uiv, où ui contient toutes 
les lettres de hauteur 1 ; et ainsi de suite, ce qui prouve que u est équivalent 
au mot obtenu en énumérant les couches de E[u). □ 

En particulier, les empilements forment un monoïde isomorphe à pc(G) ; 
notons que le produit des empilements (Ci, . . . , C„) et (C(, . . . , G'^) n'est pas 
en général égal à la suite de configurations (Ci, . . . , C^, C(, . . . , C^), qui ne 
constitue pas un empilement, sauf dans le cas exceptionnel oià C( Ç y(C„). 
Dans le langage géométrique de X.G.Viennot, on compose deux empilements 
en plaçant le second au dessus du premier, et en laissant les pièces du second 
tomber sur celles du premier. 
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Exemple 1 (suite): si on ajoute le sommet c à l'empilement E 
de la figure ^ c'est à dire si on calcule le produit E' = E.c, on 
a : 

E' = {Ci,C2,C'„C4), C!, = CsU{c} 

car c commute avec tous les éléments de C3 et C4, mais pas avec 
beC2. 

Notons aussi que lorsqu'on lit un mot u de droite à gauche, on obtient 
un nouvel empilement, dual de celui obtenu en lisant normalement le mot 
de gauche à droite ; cette transformation correspond géométriquement à une 
inversion du sens de la gravité. 

Grâce au théorème |^, les couches d'un empilement peuvent être in- 
terprétées comme une factorisation canonique d'un élément du monoïde par- 
tiellement commutatif ; c'est sous cette forme que ce théorème est apparu 
pour la première fois, dans |^, théorème 1.2], sous le nom de V-factorisation. 

3.3 Séries formelles et théorème d'inversion 

A partir du monoïde pc(G) - que nous ne distinguerons plus désormais 
du monoïde des empilements - on peut construire l'algèbre des séries formelles 
à coefficients rationnels (par exemple) sur ce monoïde : c'est l'algèbre des 
séries en variables partiellement commutatives. Un élément de cette algèbre 
est une application A qui associe à chaque empilement E (ou à sa trace) un 
coefficient A(£') dans le corps de base choisi ; on note : 

A = Y. ^(^)^ 

empilement E 

et le produit est défini naturellement par : 

^Ai(i?i)i?i I |^A2(i?2)i^2) =E ( E ^l{El)A2{E2)\ E 
El / \E2 Je \EiE2=E ) 

On vérifie facilement qu'il n'existe qu'un nombre fini de factorisations de E 
sous la forme E1E2, ce qui assure la cohérence de la définition du produit. 

Considérons en particulier les séries suivantes, qui sont en fait des polynômes 
lorsque le graphe G est fini : 

r= Y. c T= Y (-1)""^ (11) 

configuration C configuration C 

et les séries analogues (série caractéristique et série alternée) pour les em- 
pilements : 

empilement E empilement E 

en rappelant que \E\ désigne le nombre de cellules de E. 
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Exemple 2: soit G le graphe : • • • 

a h c 

Omettons les crochets [ ] pour désigner les traces des mots dans 

pc(G) ; le polynôme T vaut : 

1 — a — b — c + ac 

tandis que la série G débute par : 

1 + a + h + c + a^ + ab + ac + ha + h^ + hc + ch + c^ + ... 

(noter l'absence du terme ca, car ac = ca; un comptage manuel 
montre que la série comporte ensuite 21 termes de degré 3). 

On a le résultat fondamental suivant : 

Théorème 2 La série caractéristique des empilements, Q, et la série al- 
ternée des configurations, T, définies par ( |71| j et (jï^j, sont inverses l'une 
de l'autre. De même, G est l'inverse de T. 

Preuve: soit une configuration C et un empilement E ; l'équation CX = 
E, dans le monoïde des empilements, a une solution (unique) X si et seule- 
ment si : C Ç Cl, où Cl désigne la base de E. Donc en effectuant le produit 
r G, le coefficient d'un empilement donné E non vide vaut : E(-l)''^', la 
sommation étant étendue à tous les sous-ensembles C de Ci. Or cette somme 
est nulle (c'est la base des arguments combinatoires du type "inclusion- 
exclusion"), ce qui prouve que F G = 1. 

Le même raisonnement est valable pour l'équation : XC = E, en lisant 
les mots de droite à gauche dans pc(C) (cf. empilement dual) ce qui prouve 
que GF = 1. En changeant chaque variable en son opposée, on obtient aussi : 

rG = Gr = i. □ 

Ce théorème apparaît dans [|ï^], formulé avec des valuations ; sa source, 
citée dans |ï^, est à nouveau ||6|, théorème 2.4], oii il est présenté comme un 
calcul de fonction de Môbius (cf. la preuve) . Cas particuliers : 

- si C est le graphe complet sur A, on retrouve la formule : 

—^ = A* 
1-A 

dans l'algèbre des séries en variables non commutatives. 

- si G est le graphe totalement déconnecté - c'est à dire sans arête - sur 
A, toutes les variables commutent, la série F se factorise en : HaGAll" 
a) et la série G en : IlaeAll + o + + • • •) 
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3.4 Pyramides et théorème de décomposition 

A un empilement E est associé un ordre sur les cellules de la manière 
suivante : on dit que la cellule (a,i) précède (6, j) dans E si : 

a et 6 sont voisins, et i < j 

et on considère l'ordre engendré par cette relation. Le diagramme de Hasse 
(voir figure ^ rend cet ordre visible : c'est le graphe dont les sommets sont 
les éléments de l'ensemble partiellement ordonné considéré, et les arcs les 
couples (x, y) tels que : 

X < y et il n'existe pas d'élément z tel que x < z < y 




a c 



FiG. 4: Diagramme de Hasse de l'empilement E de la figure ^. 

Une pyramide est un empilement dont la base Ci est un singleton. Con- 
sidérons un empilement quelconque E, et une cellule 7 de ; l'ensemble des 
cellules supérieures ou égales à 7, pour l'ordre ci-dessus, forme une pyramide 
de base 7, qu'on dit engendrée par 7, et la décomposition : 

E = XP, P pyramide de base 7 (13) 

est unique. Si E est représenté par un mot u, la cellule 7 correspond à une 
occurrence d'une lettre a dans le mot u, et la pyramide P au mot formé par 
les lettres situées à droite de a dans u, et qu'aucune suite de commutations 
autorisées ne peut faire glisser à gauche de a. 

Exemple 1 (suite): l'empilement E est représenté par le mot 
u = acbegeaf, et possède les 8 décompositions suivantes (une 
par cellule), de type ([ï^), déduites de la figure ^ : 
u = cg.abeeaf = aee.cbgaf = acege.baf = acbg.eeaf 
= acbeea.gf = acbeg.eaf = acbegef.a = acbegea.f 

Aux séries caractéristiques et alternées des configurations et empile- 
ments, définies par ([ïï|) et (|I2|) , ajoutons celles des pyramides : 

n= E ^ n= E (-1)"^'^ 

pyramide P pyramide P 



15 



et les séries dérivées : 

configuration C empilement E 

(rappelons que \E\ désigne le nombre de cellules de E). 

Notons qu'il n'y a pas de pyramide vide, donc la série IT ne comporte 
pas de terme constant. 

Exemple 2 (suite): Soit G le graphe de la section p.3| ; la série 
n débute par : 

a + b + c + a^ + ab + ba + b^ + bc + cb + c^ + ... 

(suivis de 18 termes de degré 3). 

Théorème 3 La série caractéristique des pyramides est la "dérivée loga- 
rithmique droite" de celle des empilements ; plus précisément, dans l'anneau 
des séries formelles en variables partiellement commutatives, on a : 

e' = en (15) 

De même, la série alternée des pyramides est l'opposée de la "dérivée loga- 
rithmique gauche" de la série caractéristique des configurations : 

r' = -nr (16) 



Preuve: l'équation (|Î5|) résulte directement de la décomposition (p~3D, car, 
pour un empilement donné E, il existe \E\ décompositions de ce type (une 
par cellule). L'équation (|ï^ ) est une conséquence algébrique de (|^) et du 
théorème d'inversion (théorème |^), car la dérivée logarithmique de l'inverse 
de est, bien entendu, l'opposée de la dérivée logarithmique de 0. 

Précisément, en multipliant à gauche les deux termes de (15) par l'inverse 
de 0, on obtient : 

re' = n 

D'autre part, de l'égalité : 

\E1E2\E1E2 = {\Ei\ + \E2\)EiE2 = {\Ei\Ei)E2 + Ei(\E2\E2) 

on déduit par linéarité que la formule habituelle de dérivation d'un produit 
est applicable aux séries en variables partiellement commutatives, et donc : 

(re)' = r'e + rG' 

et la formule F = 1 implique donc : 

r'G = -re' d'où r'e = -n 



On en déduit (|16|) en multipliant à droite par T et en changeant les variables 
en leurs opposées, ce qui échange séries caractéristiques et alternées. □ 
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3.5 Applications aux modèles de gaz 

Corollaire 4 Soit un modèle de gaz avec particules dures, sur un graphe 
fini G, et soit : 

t = e^/'=^ 

où fi désigne le potentiel chimique du modèle, k la constante de Boltzmann, 
et T la température. Le nombre moyen de particules vaut : 

oo 
n=l 

OÙ Pn désigne le nombre de pyramides de taille n sur G. 

Preuve: Pour toute série en variables partiellement commutatives A, 
appelons projection de A la série ordinaire ô{t) obtenue en remplaçant les 
lettres de l'alphabet A par une même variable t. La projection est donc 
l'application : 

A = J2ME)E - ôit) = f: ( J2 ME)] t^ 

E n=0 \\E\=n / 

Comme par hypothèse G est fini, il n'existe qu'un nombre fini d'empilements 
de taille n, et la série ô est bien définie. 

Par définition, dériver A - cf. (14) - consiste à multiplier chaque coef- 
ficient A(E) par la taille \E\ de l'empilement ; donc la projection de A' se 
déduit de ô en remplaçant t"" par nt", autrement dit vaut tô'{t). 

Les projections de F, et H sont les séries génératrices ^{t),6{t) et Tr{t) 
qui énumèrent respectivement les configurations, empilements et pyramides 
selon leur taille ; la formule ( |Ï6|) implique : 

t7'(t) = -7r(-t)7(t) 

Or la fonction de partition Z(t) du modèle est égale à 7(t), selon la 
formule (^), section 2.3 ; et le nombre moyen de particules vaut tZ'{t)/Z{t), 



selon (|5D : le corollaire suit. □ 

Corollaire 5 Soit un modèle de gaz avec particules dures, sur un réseau G, 
et soit O un sommet de G. La densité moyenne par site vaut : 



oo 
n=l 



If-Vnt^ 



où Pn désigne le nombre de pyramides de taille n sur G, de base O. 

Preuve: sur un réseau, infini par définition, projeter une série n'a en 
général plus de sens, car il existe une infinité d'empilements de taille donnée ; 
par contre on peut projeter la série Hq des pyramides de base O, et la 
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série génératrice correspondante 7r(t) ne dépend pas du choix de O, grâce à 
l'homogénéité d'un réseau. D'autre part, soit Gk une portion finie du réseau, 
comportant k sommets (cf. section 2.4) ; on se doute que si k est suffisamment 



grand, et Gk définie de façon raisonnable, le nombre de pyramides de taille n 
sur Gfc, et de base a, dépend peu du choix de a ; autrement dit la série iTk{t) 
qui énumère toutes les pyramides sur G^, vaut approximativement kTr{t). 
On peut montrer qu'on a exactement : 

vr(t) = lim ^ (17) 

Par définition, la densité par site est égale au nombre de particules divisé 
par le nombre de sommets ; le corollaire ^ est donc une conséquence du 
précédent. □ 

Ce corollaire peut être interprété comme la limite thermodynamique du 
corollaire ^, au sens de la section 2.4 ; de même que les physiciens doivent 



définir la fonction de partition d'un gaz sur un réseau par une moyenne 
logarithmique, ici il faut remplacer l'énumération des configurations par celle 
des pyramides de base fixée, qui seule garde un sens. Notons cependant que 
notre démarche, en dehors de la formule (|ï^, est purement algébrique, et 
repose sur les deux théorèmes fondamentaux de l'algèbre des séries formelles 
en variables partiellement commutatives (inversion et décomposition). 



3.6 Empilements stricts 

Un empilement est dit strict si deux couches consécutives ne contiennent 
jamais la même lettre, autrement dit si la condition (|^) est remplacée par : 

CiÇy(Q_i)-Q_i, l<i<n 

Si u est un mot de A* qui représente l'empilement E, E est strict si et 
seulement si, pour toute décomposition : 



u = ui a U2 a U'i 



il existe une lettre de «2 qui ne commute pas avec a (voir par exemple la 
formule (|ÏÔ|) pour le prouver) ; on dira qu'un tel mot est lui aussi strict. 

Théorème 6 La série caractéristique des empilements, Q, se déduit de celle 
des empilements stricts en substituant 



a 9 
= a + a^ + a-^ + 



1 - ( 

à chaque lettre a de l'alphabet A. Inversement, on passe des empilements 
aux empilements stricts par la substitution : 



l + a 



18 



Preuve: la seconde affirmation est une conséqTicncc immédiate de la première, 
car l'inverse de la fonction a — a) est la fonction a «/(l + o.). 

Considérons ensuite un empilement E, et u un mot qui représente E ; 
écrivons u en regroupant les lettres consécutives égales : 

u = aja2^ . . . avec 7^ Oj+i pour 1 < i < m 

et parmi tous les mots u équivalents qui représentent le même empile- 
ment E, choisissons-en un avec le nombre m de facteurs minimal. Alors 
V = aia2 ■ ■ ■ ttm est strict ; sinon il existe i < j avec = aj qui commute 
avec tout Ofc, pour i < k < j, et on peut diminuer le nombre de facteurs de 
u en regroupant et aj' en a^~^^^ ■ 

Il reste à montrer que la trace du mot strict ainsi obtenu est unique ; 
soient v = 0102 . . . et v' = 6162 • • • deux mots stricts tels que u = 
a^^a2^...aj^ et u' = 6^^ 62^... 6*" soient équivalents; soit bi la première 
occurrence de ai dans v' ; puisque u ^ u' , hi commute avec pour k < i. 
Posons w = bi . . . . . . 6„ ; on a donc v' aiw. D'autre part, si par 

exemple ri > Si, il existe j > i tel que bj = bi commute avec bk pour 
i < k < j, ce qui contredit le fait que v' soit strict ; on fait un raisonnement 
symétrique si ri < s^. Donc ri = Sj, et on peut simplifier l'équivalence u ^ u' 
par o^^ ; par récurrence sur la longueur des mots, on sait que a2 ■ ■ ■ ~ w, 
d'oii V ^ v'. □ 
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4 Animaux dirigés 



4.1 Définition 

Les physiciens appellent animal un ensemble de points, ou cellules, sur 
un réseau, qui se développe à partir d'une configuration initiale par adjonc- 
tion de cellules sur les sites vides voisins de ceux occupés (deux sites sont 
voisins s'ils sont reliés par une arête du réseau). La taille d'un animal est 
égale au nombre de ses cellules. Si la configuration initiale est connexe, en 
particulier si elle est réduite à un seul point, l'animal est connexe : deux cel- 
lules distinctes peuvent toujours être reliées, dans le réseau, par un chemin 
dont tous les sommets appartiennent à l'animal. La génération de ces ani- 
maux obéit à différentes lois, selon qu'ils servent à modéliser des problèmes 
de percolation, de croissance de cristaux ou de polymères, etc. . . 

Les animaux dirigés sont contraints de se développer dans des directions 
privilégiées, définies en orientant les directions des arêtes du réseau ; pour 
chaque cellule il existe donc dans l'animal un chemin orienté depuis la con- 
figuration initiale, appelée source, vers la cellule ; on suppose en général que 
la source est réduite à un point. La figure |5| montre des exemples d'ani- 
maux dirigés sur réseaux plans ; notons que le cas (b) peut être aussi obtenu 
à partir d'un réseau carré nnn (voir figure ^ section |2.4| ), avec les direc- 
tions privilégiées Est, Nord et Nord-Est ; un tel réseau est souvent appelé 
triangulaire. 




Nous allons montrer, section 4.4, que sur certains types de réseaux, les 
animaux dirigés correspondent à des empilements. Cette équivalence passe 
par le coloriage des réseaux, qui permet de construire une variante des em- 
pilements, et de définir clairement les réseaux pour lesquels l'équivalence est 
observable ; c'est l'objet des deux sections suivantes. 
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4.2 Empilements coloriés 



Nous avons présenté, section la bijection fondamentale entre mots 
en variables partiellement commutatives, et empilements ; la clef de cette bi- 
jection réside dans la répartition en couches des lettres successives d'un mot, 
de telle sorte que l'ordre des couches spécifie entièrement l'ordre des lettres 
dans le mot, aux commutations autorisées près. Autrement dit, on ajoute 
un nouveau sommet a à un empilement en respectant la règle fondamentale 
suivante : la hauteur de a est strictement supérieure à celles de toutes les 
fibres Fh, pour b voisin de a. 

Il est possible d'introduire des variantes dans la définition des couches, 
tout en conservant cette règle ; on obtient alors de nouvelles variétés d'em- 
pilements, toutes équivalentes entre elles, puisque isomorphes au monoïde 
partiellement commutatif. La variante standard, étudiée jusqu'à présent, 
consiste à placer a aussi bas que possible - voir formule (|ï^). Dans cette sec- 
tion, nous ajouterons une contrainte, qui impose à tous les sommets d'une 
couche d'être de la même couleur. 

Un graphe colorié est un graphe G dont les sommets sont coloriés avec 
r couleurs, notées 1, 2, . . . , r, de telle sorte que deux sommets voisins soient 
de couleurs différentes ; on notera Gg l'ensemble des sommets de couleur 
•S (1 < s < r) ; la figure ^ montre un réseau linéaire et un réseau carré 
coloriés avec deux couleurs, ainsi qu'un réseau hexagonal colorié avec trois 
couleurs. 
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FiG. 6: Réseaux coloriés (a) linéaire (b) carré (c) hexagonal 

Dans un empilement colorié chaque couche Ci {1 < i < n) est un sous- 
ensemble de G s , avec i = s mod r ; les couches sont donc successivement de 
couleurs l,2,...r, et ainsi de suite de façon cyclique. La définition d'une 
coloration garantit que deux cellules d'une même couche ne sont jamais 
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voisines dans G. 

Dans un empilement standard, chaque cellule qui ne fait pas partie de 
la base, repose sur une voisine de la couche immédiatement inférieure ; dans 
un empilement colorié, on impose seulement qu'une telle cellule possède au 
moins une voisine dans l'une des r couches inférieures. Autrement dit le 
placement d'un nouveau sommet a dans un empilement colorié se fait selon 
le même principe que dans un empilement standard, mais en tenant compte 
de la couleur de a dans le choix de la couche la plus basse au-dessus des 
fibres voisines. 

Bien qu'il n'existe pas de relation simple entre les couches d'un em- 
pilement standard, et celles d'un empilement colorié, les mots obtenus par 
énumération des couches dans un empilement standard, et dans l'empilement 
colorié correspondant, sont équivalents. 

La figure |^ montre un empilement colorié sur le réseau linéaire G, et 
l'empilement standard associé ; les sommets de G sont représentés par des 
entiers : les entiers pairs sont noirs, les impairs blancs ; deux sommets i et j 
commutent si > 2. Les mots obtenus par énumération des couches sont 

respectivement 0102302302401 et 0102030203241, qui sont bien équivalents. 
La cellule la plus haute de la fibre Fi repose sur la fibre Fq dans l'empilement 
colorié, tandis qu'elle repose sur F2 dans l'empilement standard : le paradoxe 
n'est qu'apparent, car et 2 commutent. 



01234 01234 
(a) (b) 

FiG. 7: Empilements (a) colorié (b) standard, sur réseau linéaire 

4.3 Réseaux feuilletés 

A un réseau G de dimension d, colorié avec r couleurs, on peut associer 
un réseau feuilleté H, de dimension d + \, défini comme suit : les sommets 
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de H sont les couples {a,i), pour tout entier i et tout sommet a de G, de 
couleur s = i mod r ; un arc de H joint (a,i) à (6, j) si a et 6 sont voisins 
dans G, et si i<j<i + r. 

Pour une valeur donnée de l'entier i, l'ensemble des sommets (a, i) de 
H sera appelé couche numéro i ; par définition, une couche est monocolore. 
Exemples : 

1. Le réseau feuilleté associé au réseau linéaire colorié de la figure |^ (a), 
est le réseau carré de la figure]^ (a) (les couches sont horizontales). 






FiG. 8: Réseaux feuilletés 

2. Le réseau feuilleté associé au réseau carré colorié de la figure |6|-(b), est 
le réseau cubique centré ; on l'obtient en ajoutant aux sommets d'un 
réseau cubique, coloriés en blanc, les centres des cubes élémentaires, 
coloriés en noir : voir la figure |^(b). Ainsi chaque sommet de la couche i 
possède 8 voisins immédiats (les sommets du cube dont il est le centre), 
répartis dans les couches i — 1 et i + 1. 

3. Le réseau feuilleté associé au réseau hexagonal colorié de la figure |6|- 
(c), est le réseau cubique nnn - voir section 2.4 pour la définition des 
réseaux nnn. Chaque hexagone du réseau plan est considéré comme 
la projection d'un cube "posé sur la pointe" ; chaque sommet de la 
couche i, possède 3 voisins immédiats dans la couche i + et 3 "voisins 
suivants" dans la couche i + 2 : voir la figure ^-(c). 

Un réseau feuilleté peut être orienté en dirigeant les arcs des couches 
inférieures vers les couches supérieures (il n'y a jamais d'arête à l'intérieur 
d'une couche) ; il peut être étendu en ajoutant les arcs (a, i) — > (a, i + r). Par 
exemple le réseau feuilleté étendu associé au réseau linéaire est un réseau 
triangulaire, obtenu en ajoutant des arcs verticaux à la figure ^(a). 



4.4 Equivalence entre animaux dirigés et empilements 

Théorème 7 Soit G un réseau colorié, de dimension d, H le réseau feuilleté 
orienté issu de G, de dimension d+1, et H' le réseau étendu. Il existe une 
bijection entre pyramides (resp. pyramides strictes) construites sur G, et 
animaux dirigés, de sources ponctuelles, sur H' (resp. H). 
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Preuve: soit r le nombre de couleurs de G ; par définition d'un empilement 
colorié E, chaque cellule {b,j) qui n'appartient pas à la base de l'empilement, 
possède une voisine (a, i) dans l'une des r couches inférieures. Si j = i + r, 
alors a = b, car les couches i et i + r ont même couleur ; donc si l'empilement 
est strict, on a j < i + r. Par définition d'un réseau feuilleté orienté, il existe 
un arc {a,i) {b,j) dans H', et même dans H si l'empilement est strict. 

De proche en proche, on en déduit que, pour toute cellule de E, il existe 
un chemin dans H (ou dans H') depuis la base de E jusqu'à cette cellule. 
Un empilement colorié sur G est donc un animal dirigé sur H', et même sur 
H s'il s'agit d'un empilement strict ; et vice-versa. Si la source de l'animal 
est ponctuelle, la base de l'empilement est réduite à un point : c'est donc 
une pyramide. 



Enfin on a vu section 4.2 qu'il existe une bijection naturelle entre em- 



pilements standard et empilements coloriés. □ 
Exemples : 

1. Lorsque le nombre de couleurs r vaut 2, il n'existe aucune différence en- 
tre empilements stricts standard et empilements stricts coloriés. Donc 
l'animal dirigé de la figure |l-(a) est exactement une pyramide stricte 
sur réseau linéaire : il suffit de faire subir au dessin une rotation de 
7r/4 dans le sens trigonométrique pour le vérifier. 

2. L'animal dirigé de la figure |^(b) est une pyramide coloriée (non stricte) 
sur réseau linéaire : il faut faire subir cette fois au dessin une rotation de 
7r/6 dans le sens trigonométrique. La pyramide coloriée de la figure]^ 
(a) est un autre exemple d'animal dirigé sur réseau triangulaire ; la 
pyramide standard de la figure 0-(b) est donc la pyramide associée à 
cet animal. 

3. Une pyramide stricte sur réseau carré correspond à un animal dirigé 
sur un réseau cubique centré ; l'animal peut se développer, dans un 
repère convenablement choisi, selon l'un des 4 vecteurs : 

(1,1,1), (1,-1,1), (-1,1,1), (-1,-1,1) 

auxquels il faut ajouter le vecteur (0, 0, 2) pour traiter les pyramides 
quelconques ; voir figure ^(b). 

4. Une pyramide stricte sur réseau hexagonal correspond à un animal 
dirigé sur un réseau cubique nnn ; l'animal peut se développer selon 
l'un des 6 vecteurs : 

(1,0,0), (0,1,0), (0,0,1), (1,1,0), (1,0,1), (0,1,1) 

auxquels il faut ajouter le vecteur (1, 1, 1) pour traiter les pyramides 
quelconques. Noter qu'avec le repère choisi ici, les couches ne sont plus 
horizontales, comme sur la figure ^-(c), mais d'équations x + y + z = i. 
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En combinant ce théorème, dû à X.G.Viennot, avec le corollaire |5| de 



la section 3.5, et avec le théorème g, qui relie les séries caractéristiques des 



empilements selon qu'ils sont stricts ou non, on retrouve le résultat suivant 



de D.Dhar llO|] 



Corollaire 8 Le problème de Vénumération des animaux dirigés sur un 
réseau feuilleté H issu d'un réseau G, est équivalent à celui du calcul de 
la densité moyenne d'un gaz à particules dures sur G, en fonction de la 
température. □ 

Il est intéressant de comparer les méthodes employées en physique théo- 
rique et celles, combinatoires, présentées dans cet article. Dans [|ÏÔ|, D.Dhar 
emploie l'approche suivante : 

" Un site j est appelé successeur du site i s'il existe un arc 
de i vers j. Pour un ensemble de sites G, l'ensemble des suc- 
cesseurs sera noté S (G). (...) Définissons pour chaque site i du 
réseau un temps t{i), à valeurs entières, tel que t{j) > t{i) + 
1 chaque fois que j est un successeur de i. (Dans un animal 
dirigé) les configurations autorisées de sites occupés, sur la sur- 
face t = r, ne dépendent que des configurations aux instants 
précédents t < T — 1. Considérons d'abord, pour simplifier, le 
cas oir elles dépendent seulement de la configuration sur la sur- 
face t = T — 1 (cas du réseau carré, et du réseau cubique centré). 
Soit G l'ensemble des sites occupés sur une hyper-surface t = c^^^. 
Définissons les fonctions génératrices Aq (x) , qui énumèrent selon 
leur taille les animaux dirigés de source G. 

La propriété markovienne des animaux dirigés implique que ces 
fonctions génératrices vérifient les récurrences : 



où la sommation est étendue à tous les sous-ensembles non vides 
D de S (G). Notons que le temps sur les sites de D est supérieur 
de 1 au temps sur les sites de C. Il est facile de généraliser ces 
récurrences au cas où les configurations autorisées à l'instant t = 
T dépendent aussi des configurations aux instants t <T — 2 (cas 
des réseaux nnn, etc. . . ). Ces équations forment la hiérarchie de 
Bogoliubov-Born-Green-Kirkwood-Yvon (sic) pour le problème 
de rénumération des animaux dirigés. " 

D.Dhar compare alors ces équations à celles qu'on obtient en simulant 
l'évolution d'un gaz : à chaque incrémentation du temps, chaque site évolue 
aléatoirement, la probabilité d'occupation d'un site donné ne dépendant que 
de la configuration des sites voisins, avec une formule simple déduite de la loi 
de Gibbs dsh ; lorsque i — > oo, le gaz approche ainsi de l'équilibre thermique. 
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Ce principe de simulation, classique en physique statistique, est aussi à la 
base des réseaux de neurones (voir par exemple [p7||). 

" (Etudions) l'évolution stochastique dans le temps d'un gaz 
sur un réseau de dimension d — 1. Ce réseau est constitué de r 
sous-réseaux entrelacés Li, L2, ■ ■ ■ , Lr, chacun d'entre eux étant 
isomorphe à un hyperplan r = c^^^ du réseau de dimension d 
(sur lequel se développent les animaux dirigés). (...) A l'instant 
r = mr + i (1 < i < r), la configuration des sites occupés sur le 
sous-réseau Li évolue aléatoirement, tandis que les autres sous- 
réseaux restent inchangés. Lorsque r — > 00, la probabilité des 
différentes configurations du gaz tend vers une distribution lim- 
ite invariante. Pour des lois de transition arbitraires, le calcul de 
la distribution invariante est tout à fait difficile. (Ici) les lois de 
transition sont particulièrement simples, et correspondent aux 
algorithmes de Monte-Carlo utilisés couramment pour étudier 
les lois d'équilibre d'un gaz sur un réseau avec particules dures 
(...) La probabilité d'occupation d'un site i vaut si l'un des 
sites successeurs de i est occupé ; et elle vaut p s'ils sont tous 
vides. " 

D.Dhar remarque qu'il suffit de "renverser le temps", c'est à dire de poser 
r = —t, pour faire disparaître le paradoxe qui fait dépendre la probabilité 
d'occupation d'un site i de la configuration de ses successeurs, au lieu de ses 
prédécesseurs. Il écrit : 

" Soit P{C) la probabilité que tous les sites de l'ensemble C, 
inclus dans une hyper-surface r = c^^^, soient occupés. Je suppose 
r suffisamment grand pour que P{C) ne dépende pas de r. Cette 
probabilité est égale à p'*^' multiplié par la probabilité que tous 
les sites de S{C) soient vides. En utilisant le principe d'inclusion- 
exclusion 

P(C) =pl^l |^l + ^(-l)l^lp(D)^ 

oii la sommation est étendue à tous les sous-ensembles non vides 
D de S{C). On voit donc que 

Ac{x = -p) = {-lf\p{C). 

Ainsi la fonction génératrice des animaux dirigés issus d'un seul 
point source (sur un réseau de dimension d) est déterminée par 
la densité moyenne de sites occupés dans un gaz à l'équilibre sur 
un réseau de dimension d — 1. " 

On voit comment le théorème d'inversion est remplacé par l'utilisation 
du principe d'inclusion-exclusion en calcul des probabilités. L'introduction 
de sous-réseaux entrelacés, qui a été utile aussi dans le présent article pour 
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éclaicir la bijection entre empilements et animaux, n'est pas expliquée par 
D.Dhar. Elle est rendue nécessaire par le modèle des particules dures, car 
si à l'instant r tous les sites changent aléatoirement d'état, on peut voir 
apparaître des particules sur des sites voisins, chacun étant pourtant entouré 
de sites vides à l'instant précédent. 

Nos méthodes algébriques sont donc parentes de celles des physiciens, 
dont on ne peut qu'admirer l'ingéniosité ; mais l'algèbre et la combinatoire, 
substituées aux récurrences déduites du calcul des probabilités, permettent 
d'éliminer l'imprécision de certains arguments ; surtout elles permettent une 
étude plus fine des objets : dans la section suivante, l'algèbre des empilements 
nous donnera des décompositions canoniques des animaux en dimension 2, 
qui permettent à la fois de retrouver très simplement les séries génératrices 
qui les énumèrent, et de construire des animaux aléatoires, avec distribution 
uniforme, en temps linéaire. 
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5 Solution du modèle linéaire 



5.1 Décomposition en équerres 



Le cas le plus simple de modèle de gaz à particules dures sur un réseau 
G, est celui où G est de dimension 1, autrement dit où G est une chaîne 



infinie, isomorphe à Z. D'après la section 4.4, les pyramides strictes sur G 
sont exactement les animaux dirigés, de source ponctuelle, sur réseau carré. 

Repérons les sommets de G par des entiers, et considérons une pyramide 
P de base 0. On peut décomposer P en remarquant que : 

- si P ne contient aucune cellule au-dessus de 1, autrement dit si la fibre 
Fi est vide, alors toutes les fibres Fi, pour i > 1, sont vides aussi; 
on dira dans ce cas que P est une équerre de base 0. Cette propriété 
repose sur une particularité évidente du modèle linéaire : il n'existe 
pas de chemin de vers i > 1, qui ne passe pas par 1. 

- si Fi n'est pas vide, soit (1, k) la cellule la plus basse sur cette fibre, et 
soit Q la pyramide engendrée par cette cellule ; la décomposition (lïî 



de la section 3.4 devient : 



(18) 



P = LQ 

où L est une équerre, car Q contient la fibre Fi. 
La décomposition ( p^ ) peut évidemment être itérée, en décomposant Q, 
etc. . . ; elle est schématisée, ainsi que son itération, sur la figure La plus 
grande valeur de i pour laquelle la fibre Fi n'est pas vide, est appelée demi- 
largeur droite de la pyramide (ou de l'animal) P ; le nombre d'équerres qui 
interviennent dans la décomposition de P est égal h i + 1; une équerre est 
une pyramide de demi-largeur droite nulle. 





••■■-1012 -1 1 •■ i ■■■■ 

FiG. 9: Schéma de décomposition d'une pyramide en produit d'équerres 

La figure |ÏÔ| illustre ces décompositions sur l'exemple de l'animal de la 
figure |5|-(a), dont la demi-largeur droite vaut 4; la partie (a) de la figure 
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montre la première décomposition, avec les cellules de L marquées en noir, 
et celles de Q en blanc ; la partie (b) montre la décomposition itérée de 
Q comme produit de quatre équerres. On notera que la pyramide Q ne 
semble pas former un animal sur la figure ^^(a) ; (b) montre Q sous forme 
d'animal, après la "chute" de quatre cellules, une fois l'équerre L disparue. 
Cet exemple montre combien l'algèbre des empilements est essentielle pour 
ces décompositions, invisibles si l'on s'en tient à la définition brute d'un 
animal. 




(a) (b) 
FiG. 10: Exemple de décomposition d'une pyramide en produit d'équerres 

En notant Ilj (resp. Aj) la série caractéristique des pyramides (resp. 



équerres) de base i, la décomposition (|18D devient (en tenant compte du cas 
où P est déjà une équerre) : 

no = Ao(i + ni) (19) 

Une équerre L de taille strictement supérieure à 1, peut être à son tour 
décomposée ; soit x la base de L : 

- si X est seule sur la fibre Fq, alors la pyramide M, engendrée par la 
cellule (—1, 1), est une équerre, et : 

L = xM (20) 

- sinon, soit (0, /c) la cellule la plus basse sur Fq, avec k > 1; cette 
cellule engendre une pyramide A^, qui est à son tour une équerre, et la 
décomposition : 

L = xMN (21) 

laisse une pyramide M, de base (—1,1), dont on a vidé la fibre Fq ; M 
est donc une équerre. 
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Ces décompositions sont schématisées sur la figure 11, parties (a) et (b) ; 
la partie (c) illustre (b) sur l'exemple de l'équerre noire de la figure |ï^(a) ; 
l'équerre M est laissée noire, et les cellules de apparaissent en blanc. On 
peut voir à nouveau que cette décomposition n'est pas visible sans considérer 
un animal comme un empilement, car n'est pas connexe. 






• ■ • • -1 1 

(a) (b) (c) 

FiG. 11: Décomposition d'une équerre en produit d'équerres 

Les décompositions ( pÔ|) et (^) entraînent, pour les séries caractéristi- 
ques, en tenant compte du cas où L = x (pyramide réduite à sa base) : 



An 



x(l + A_i + A„i Ao 



(22) 



Les équations (pï^ ) et ( ^4 ) concernent les pyramides strictes, ou les an- 
imaux dirigés sur réseau carré ; pour les pyramides générales, qui corre- 
spondent aux animaux dirigés sur réseau triangulaire, on peut appliquer le 
théorème ^ ou remarquer que ( |T9| ) reste valable, et que (22) devient : 



Ao = x(l + A_i + Ao + A_i Aq) 



(23) 



car il y a une seconde décomposition de type |TÏ|-(a), avec M, de base 0, 
reposant verticalement sur x. 



5.2 Séries génératrices 

On obtient les séries génératrices des pyramides et des équerres en proje- 
tant (voir section ^) les équations obtenues pour les séries caractéristiques 
dans la section précédente ; ces projections ne dépendent pas de la base. 
En notant iï{t) (resp. \{t)) la série génératrice qui énumère les pyramides 
strictes (resp. les équerres strictes), (|Ï9|) et (|^) fournissent : 

7r(i) = A(t)[l + 7r(t)] (24) 
A(t) = t[l + \{t) + \{tf] (25) 
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Ces équations sont facilement résolues : 



A(t) = (26) 

= t + f + 2t^ + At^ + 9t^ + 21t^ + blf + 127t^ + 323t^ + . . . 




-it) = ^h/l^-l| (27) 
= t + 2f2 + 5t^ + 13t^ + 35t^ + 96f + 267 f + 750t^ + 2123t^ + . . . 
Pour énumérer les équerres et les pyramides générales, il faut remplacer 



(25) par l'équation déduite de (|25 

X{t) =t[l + 2X{t) + X{tf] (28) 

et on obtient cette fois : 



= t + 2*2 + 5t^ + Ut^ + A2t^ + 132t^ + 429*^^ + 1430t® + . . . 



= t + 3*2 + lOt^ + 35t^ + 126t^ + 462t^ + 1716t'^ + 6435t® + . . . 



Il n'existe pas de forme close pour les coefficients des séries ( Pq) et (27), 
contrairement au cas des séries (|29| ) et (^Ôj). On vérifie que, conformément 
au théorème ^, on peut passer du premier couple de séries au second par 
la substitution t ^ t/{l — t), et inversement du second au premier par la 
substitution t ^ t/{l + t). 

En appliquant le théorème ^, section 4.4, on retrouve les résultats de 
|Ô|, [H, II, 0] : 



Théorème 9 La série ^21{ ) - resp. I[3(\ ) - est la série énumératrice des 
animaux dirigés, de source ponctuelle, sur réseau carré - resp. triangulaire - 
, selon leur taille. □ 



Si, d'autre part, on applique le corollaire |5| de la section 3^, on obtient : 



Théorème 10 La densité moyenne d'un gaz avec particules dures, sur ré- 
seau linéaire, vaut : 

d{t) = -(l- , ^ 
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Preuve: d'après le corollaire ^, cette densité vaut — 7r(— t), où 7r(t) est la 
série énumératrice des pyramides, donnée par (^Ô|). □ 

Rappelons que la variable t est liée à la température absolue T par 
t = e'^^^^ , où /X désigne le potentiel chimique, et k la constante de Boltzmann 
(cf. section p.3p . Si le potentiel chimique est positif, l'énergie diminue lorsque 
le nombre de particules augmente ; lorsque la température T croît de à +00, 
t décroît de +00 à 1, et d{t) décroît de 1/2 (configuration d'énergie minimale : 
un site sur deux est occupé) à (1 — l/\/5)/2 (toutes les configurations sont 
équiprobables). On note que la densité ne présente pas de singularité, bien 
que son développement en série ait pour rayon de convergence 1/4. 

D.Dhar, dans |ï^, déduit le théorème ^ du théorème |ï^, bien connu des 
physiciens, alors que notre approche donne les deux résultats simultanément. 



5.3 Chemins de Dyck et chemins de Motzkin 

Les formules de la section précédente sont classiques en combinatoire, 
dans d'autres contextes. Rappelons qu'un chemin de Dyck est un chemin 
dans le plan discret, composé de pas ascendants (1,1) et de pas descendants 
(1,-1), et qui joint l'origine (0,0) à un point de l'axe des abscisses, sans 
quitter le quart de plan supérieur : voir figure |ï^(a). Un chemin de Motzkin 
est défini de façon semblable, mais comporte en outre des pas horizontaux 
(1,0) : voir figure ^{h). 

Dans un chemin de Motzkin r-colorié, il y a r sortes de pas horizontaux ; 
dans la suite r = 0, 1 ou 2 : les chemins 0-coloriés sont les chemins de Dyck, 
les chemins 1-coloriés sont les chemins de Motzkin ordinaires, et les chemins 
2-coloriés seront dits hicoloriés. 

Un préfixe de Motzkin est un chemin de Motzkin éventuellement incom- 
plet, c'est à dire qui ne se termine pas forcément sur l'axe des abscisses ; la 
hauteur d'un préfixe désigne l'ordonnée du point d'arrivée. 
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(a) (b) 
FiG. 12: Chemins (a) de Dyck (b) de Motzkin 



La longueur d'un chemin compte le nombre de pas ; un chemin peut être 
vide : il est alors de longueur nulle. En notant un pas ascendant a, un pas 
descendant 6, et un pas horizontal de couleur i , Cj (1 < i < r), on associe à 
un chemin de Motzkin r-colorié un mot de Motzkin, sur un alphabet kr + 2 
lettres ; les mots de Dyck correspondent au cas r = ; le pas ci sera parfois 
noté simplement c. 
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Un chemin de Motzkin non vide W se décompose selon le schéma de la 
figure |ï^, en W = CiU - cas (a) - ou en = allbV - cas (b) ; U et V 
désignent à nouveau des chemins de Motzkin (éventuellement vides) : on se 
trouve dans le cas (a) ou (b) selon que W débute par un pas horizontal ou 
ascendant ; dans ce dernier cas, V débute dès le premier retour sur l'axe des 
abscisses, d'oii l'unicité de la décomposition. 





(a) (b) 

FiG. 13: Décompositions de chemins 



En désignant par M le langage des mots de Motzkin r-coloriés, on en 
déduit l'équation : 

M = 1 + ^ CiM + aMbM (31) 



l<î<r 

(oii la somme J2 disparaît pour r = 0). La série énumératrice des mots 
de Motzkin selon leur longueur vérifie donc : 

/x(t) = 1 +rtfi{t)+ t'unit f 

dont les solutions sont ( p6|) et ( p9| ) pour r = 1 et r = 2, avec X{t) = tfi{t). 
Un préfixe W d'un chemin de Motzkin, ou bien est un chemin de Motzkin 



complet, ou bien se décompose selon la figure 14, en un chemin de Motzkin 
U suivi d'un préfixe V, soit W = UaV. Cette décomposition est unique, 
car U se termine au dernier passage sur l'axe des abscisses. Donc la série 
énumératrice p(t) des préfixes est liée à celle des chemins complets par : 



pit)=fi{t)[l + tp{t)] 




FiG. 14: Décomposition d'un préfixe 

Nous appellerons célibataire un pas ascendant conduisant du niveau y au 
niveau y + 1 - ie. un pas (x, y) — > (x + 1, y + 1) -, et qui n'est suivi d'aucun 
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pas descendant ramenant le chemin au niveau y. Si on itère la décomposition 
ci-dessus d'un préfixe W, on obtient : 



W = UoaUia...Ui (32) 

où les occurrences de a désignent les pas ascendants célibataires, et les Ui des 
mots de Motzkin (éventuellement vides). Cette décomposition est unique, 
et la hauteur du chemin est égale au nombre l de pas célibataires. 

L'analogie entre les décompositions de chemins et d'animaux, fournit 
une bijection résumée dans le théorème suivant, dû à (mais avec une 
construction totalement difi'érente) : 

Théorème 11 // existe une bijection entre, d'une part : animaux dirigés, de 
source ponctuelle, sur réseau carré (resp. triangulaire), de taille n, et d'autre 
part : préfixes de Motzkin (resp. préfixes bicoloriés), de longueur n — 1. 
Dans cette bijection, les animaux de demi-largeur droite k correspondent 
aux chemins de hauteur k ; en particulier les équerres correspondent aux 
chemins de Motzkin. 

Preuve: la bijection se construit de manière récursive ; détaillons-la dans 
le sens chemin — > animal, et notons-la f3 ; rappelons que |Î7| désigne la 
longueur d'un chemin U, et |P| la taille d'un animal P. Au chemin vide 
(de longueur nulle) on associe l'animal réduit à une cellule (de taille 1). 

Soit W un chemin de Motzkin non vide; si W = CiU, on associe récur- 
sivement kU un animal en équerre M = f3{U), et on pose P(W) = L = xM, 
avec les notations de ( [2Ô| ) ; si W est un chemin de Motzkin ordinaire, la 
base de M est translatée horizontalement de -1 par rapport à x, base de L, 



comme sur la figure 11 ; si VF est un chemin de Motzkin bicolorié, la base de 
M est translatée de -1 ou de selon que W débute par ci ou C2 : on obtient 
alors un animal sur réseau triangulaire. 

De même, siW = aUbV , on associe récursivement kU qïV des équerres 
M = f3{U) et N = I3{V), et on pose (3{W) = L = xMN, avec les notations 
de (]2Ï]). La base de M est translatée horizontalement de -1 par rapport à x, 
base de L, et A'' a sa base sur la même fibre que x. 

Si W est un préfixe de Motzkin qui ne retourne pas sur l'axe des abscisses, 
alors W = U aV, où a désigne un pas ascendant célibataire ; on associe 
récursivement k U et V une équerre L = P{U) et un animal Q = f3{V), 
et on pose P{W) = P = LQ, avec les notations de (|ï^). La base de Q est 
translatée horizontalement de +1 par rappport à celle de L. 

On vérifie que ces constructions conservent la relation |/3(1^)| = \W\ + 
1 et que la demi-largeur droite de l3{W) s'accroît en même temps que la 
hauteur de W (rappelons que la hauteur de W est définie comme l'ordonnée 
du point d'arrivée). □ 



Cette bijection est écrite en langage C dans la section S.l, et la section 3.2 



contient des exemples (avec ci noté c, et C2 noté d) . Le théorème 1 1 confirme 
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les valeurs trouvées pour les coefficients des séries ( p9| ) et (pÔ|), grâce à la 
proposition suivante : 

Proposition 12 Le nombre de mots de Motzkin bicoloriés de longueur n — 1 
est égal au nombre de mots de Dyck de longueur 2n, donné par le nombre 
de Catalan : 

C _ 1 /^2n\_ 1 (2n+l\ 
" n + l\ n j 2n + l\ n j ' 

Le nombre de préfixes de Motzkin bicoloriés de longueur n — 1 est égal au 
nombre de préfixes de Dyck de longueur 2n — 1, donné par le coefficient 
binomial : 

2n-l\ _ 1 / 2n 
n j 2 \ n 

Preuve: cette proposition est classique en combinatoire ; la première par- 
tie découle simplement de l'identité entre séries génératrices : 

[l + \{t)f = l + 2\{t) + \{tf . 

Nous en donnons une preuve bijective pour être complets. Pour la pre- 
mière affirmation, il suffit d'associer à un mot de Motzkin bicolorié W , un 
mot de Dyck non vide W défini de façon récursive - e désigne le chemin 
vide : 

e^ab , ciU ^ abU' , C2U aU'b , aUbV aU'bV . 

C'est une bijection, car l'équation D = 1 + aDbD, vérifiée par la série car- 
actéristique des mots de Dyck - voir (^) -, peut s'écrire : 

D-l = ab + ab{D - 1) + a{D - 1)6 + a{D - l)b{D - 1) . 

Pour la seconde partie de la proposition, on part deW = U^aUia . . .Ui- 
voir ( |32|) -, on associe à chaque mot de Motzkin bicolorié Ui un mot de Dyck 
Ul comme ci-dessus, on appelle U'^' le mot U[ privé de son pas descendant 
final 6, et on associe à W le mot : 

W" = UliaU'{a...U[' . 

C'est une bijection car on peut retrouver la décomposition ci-dessus de W" 
en repérant ses pas ascendants célibataires de rang pair. On vérifie aussi que 
\W"\ = 2\W\ + 1. 

Enfin, un préfixe de Dyck de longueur impaire 2n — 1 possède un nombre 
impair 2A; — 1 de pas ascendants célibataires ; en remplaçant les k premiers 
par des pas descendants, on obtient un chemin joignant l'origine au point 
(2n — 1, —1). C'est encore une bijection, car les pas descendants ainsi intro- 
duits sont les pas descendants célibataires : un pas descendant conduisant 
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du niveau y + 1 au niveau y - ie. un pas (a;, y + 1) (x + 1, y) -, est appelé 
célibataire s'il n'est précédé d'aucun pas ascendant conduisant du niveau y 
au niveau y + 1, autrement dit s'il est le premier atteignant le niveau y. Ceci 
explique que les préfixes de Dyck soient comptés par un coefficient binomial. 
□ 

La technique qui consiste à repérer, dans un chemin (ou mot) quelconque 
W , les pas célibataires, conduit à la factorisation de Catalan : 

W = UobUib...UkaUk+ia...Uk+i (33) 

qui comporte k (resp. /) pas célibataires descendants (resp. ascendants), et 
est unique : voir figure |ï^. 

y 



Uo 



k+r 



b\ 



FiG. 15: Factorisation de Catalan d'un chemin quelconque 

Les préfixes de Motzkin sont les chemins sans pas célibataire descendant 
(/c = 0), et le nombre l de pas célibataires ascendants est égal à la hauteur 
du chemin ; les chemins de Motzkin sont ceux qui ne possèdent aucun pas 
célibataire, ascendant ou descendant (A; = / = 0). 

A un préfixe de Motzkin de hauteur /, on peut associer l + 1 mots en 
remplaçant les j (0 < j < l) premiers pas célibataires ascendants par des pas 
descendants; on en déduit, comme dans [11|, que la moyenne des hauteurs 
des préfixes de Motzkin r-coloriés de longueur n vaut : 

(r + 2)" 



1 



où rrin désigne le nombre de ces préfixes ; comme la largeur moyenne d'un 
animal est évidemment le double de la demi-largeur droite moyenne, le 



théorème 11 fournit le corollaire suivant, emprunté à |l]l 



Corollaire 13 Soit an le nombre d'animaux de source ponctuelle, sur réseau 
carré ou triangulaire, de taille n ( donné par le théorème ^ ; la largeur 
moyenne de ces animaux vaut : 

(r + 2)"-i 
2 ^ ' 2 

an 

avec r = 1 sur réseau carré, r = 2 sur réseau triangulaire. □ 
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Appelons compacte une configuration, sur un réseau, aussi dense que 
possible (en respectant bien entendu la règle d'exclusion entre sites voisins) ; 
sur un réseau linéaire, la distance entre deux cellules d'une configuration 
compacte est donc égale à deux. Un animal de source compacte correspond 
à un empilement de base compacte. La factorisation de Catalan fournit 
immédiatement une bijection entre mots et animaux de source compacte : 
à chaque facteur Ui de (^), on associe par le théorème 11 un animal en 
équerre Lj, et à VF on associe le produit des équerres, en espaçant de deux 
unités les bases de Lq . . . Lk, puis seulement d'une unité les bases de . . . L/. 
Voir une illustration section 6.2. Nous obtenons ainsi une nouvelle preuve 
de l'étonnant théorème découvert par [^] : 

Corollaire 14 Le nombre d'animaux de taille n et de source compacte sur 
réseau carré (resp. triangulaire) est égal à 3"""^ (resp. □ 
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6 Animaux plans aléatoires 



6.1 Code 



Les théorèmes 11 et 14 de la section 5.3 fournissent des bijections ex- 
plicites entre animaux et mots sur un alphabet à r + 2 lettres (r = 1 pour un 
réseau carré, r = 2 pour un réseau triangulaire) ; ils permettent de construire 
facilement un programme de génération aléatoire d'animaux plans, avec dis- 
tribution uniforme. La complexité de l'algorithme est linéaire en moyenne : 
le temps d'exécution est proportionnel, en moyenne, à la taille n de l'animal. 
La mémoire utilisée est proportionnelle à n, et on ne peut évidemment faire 
mieux si on mémorise l'animal ; sinon on pourrait se contenter d'une taille 
mémoire proportionnelle, en moyenne, à ^/n. Nous donnons explicitement 
ce programme, sous forme de fonctions écrites en langage C. 

static int nb_tirages; char Mot [TAILLE_MAX] ; 
static char hasard (n) 

int n; {return (char) ( 'a' + lrand48() 7» n );} 



static void mot (n, r) 
int n, r; {int i; 

for (i = 0; i < n; i++) Mot [i] = hasard (r + 2) ; 

Mot [n] = '\0'; 

} 

static void préfixe (n, r) 
int n, r; {int i, h; char pas; 

for (nb_tirages =i=h=0; i<n; nb_tirages++) { 
Mot [i++] = pas = hasard (r + 2) ; 
if (pas == 'a') h++; 
else if (pas == 'b') 

if ( — h <0)i=h=0; /*on recommence tout */ 

> 

Mot [n] = '\0'; 

> 



FiG. 16: Code pour la génération aléatoire de préfixes de Motzkin 

Il faut commencer par construire un mot au hasard, de longueur n, sur 
un alphabet de r + 2 lettres. Si l'on garde n'importe quel mot (fourni par 
une suite de tirages de lettres, effectués à l'aide d'un générateur de nombres 
aléatoires de la bibliothèque numérique Unix standard), l'animal correspon- 
dant a une source compacte (voir théorème [î^ ) . Si l'on souhaite construire 
un animal de source ponctuelle, il ne faut garder que les préfixes de Motzkin 



(voir théorème 11) ; dans ce cas on recommence la procédure depuis le début, 



dès que le chemin associé traverse l'axe des abscisses ; cet algorithme très 
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simple ne détruit pas le caractère uniforme de la distribution, et il est montré 
dans Q que sa complexité est en moyenne 2n - voir aussi 

La fonction préfixe (n, r), qui calcule dans la variable globale Mot, un 
préfixe de Motzkin aléatoire de longueur n, sur un alphabet à r + 2 lettres, 
est donnée figure La variable globale nb_tirages permet de calculer la 
complexité (expérimentale) de l'algorithme, et de vérifier, si on le souhaite, 
les prévisions de ||2|. Nous donnons aussi la fonction mot (n, r), qui fournit 
un mot aléatoire quelconque, pour être complets. Le générateur de nombres 
aléatoires utilisé, Irand^S (), doit en principe être initialisé avant emploi, par 
appel à srand48 () : voir le manuel Unix. 

static celibataires_ascendants (n) 
int n; {int i, h = 0, hmin = 0; char pas; 
for (i = n-1; i >= 0; i — ) { 
pas = Mot [i] ; 
if (pas == 'b') h++; 
else if (pas == 'a') 

if (—h < hmin) { hmin = h; Mot [i] = 'A'; } 

} 

> 

static celibataires_descendants (n) 
int n; {int i, h = 0, hmin = 0; char pas; 
for (i = 0; i < n; i++) { 
pas = Mot [i] ; 
if (pas == 'a') h++; 
else if (pas == 'b') 

if (—h < hmin) { hmin = h; Mot [i] = 'B'; } 

} 

} 



FiG. 17: Code pour marquer les célibataires 

La construction de la bijection entre mots et animaux, fait jouer un rôle 
tout à fait particulier aux pas célibataires (voir définitions pages ^ et |35| ) . 
Ici on marque les pas célibataires en transformant la lettre correspondante 
en majuscule. Nous incluons ces deux fonctions élémentaires, pour être ex- 



plicites, figure 17. 



La fonction equerre (i), sur la figure 18, construit un animal en équerre 
dont la base a pour abscisse i, autrement dit est située sur la fibre i. Les 
hauteurs des fibres sont conservées dans la variable fibre, et la fonction em- 
piler (i), détaillée plus loin, ajoute une cellule sur la fibre numéro i de l'ani- 
mal. Le reste du code de la fonction equerre est une traduction directe de la 
preuve du théorème ; les pas c et d sont les deux types de pas horizontaux, 
notés Cl et C2 section le pas d n'est utilisé que sur réseau triangulaire. On 



notera l'absence de référence explicite au pas 6, qui sert de séparateur entre 
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static int 
compteur, 

fibre [FIBRE_MAX] ; 



/* compte les lettres */ 
/* hauteurs des fibres */ 



static void equerre (i) 
int i; { 

empiler (i) ; 

switch (Mot [compteur++] ) { 
case 'a' : 

equerre (i - 1); equerre (i) ; break; 
case ' c ' : 

equerre (i - 1) ; break; 
case ' d' : 

equerre (i) ; break; 

} 

} 

void animal (n, triangulaire, compact) 
int n, triangulaire, compact; -[int i, r; 

if (triangulaire) r = 2; else r = 1; 

for (i = 0; i < FIBRE_MAX; i++) fibre [i] = -1; 

if (compact) { 

mot (n - 1 , r) ; 

celibataires_descendants (n - 1) ; 

} 

else 

préfixe (n - 1 , r) ; 
celibataires_ascendants (n - 1) ; 

for (i = compteur = 0; compteur < n; i++) { 

equerre (i) ; if (Mot [compteur - 1] == 'B') i++; 

} 

} 

FiG. 18: Code pour transformer un mot en animal 
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deux équerres : après lecture d'un pas ascendant a, et construction récursive 
d'une première équerre, par l'appel equerre on trouve forcément un 

pas descendant b, sinon le pas ascendant serait célibataire, et noté A ; on 
appelle alors equerre (i). 

La fonction animal (n, triangulaire, compact) est la fonction principale, 
qui construit un animal de taille n, sur réseau carré ou triangulaire, selon 
la valeur du paramètre booléen triangulaire ; la source est ponctuelle ou 
compacte selon la valeur du paramètre booléen compact. Dans cette fonction, 
on initialise la taille de l'alphabet selon la nature du réseau, on rase les fibres, 
et on construit un mot aléatoire quelconque, ou un préfixe de Motzkin, selon 
que la source de l'animal est compacte ou ponctuelle. Après transformation 
des pas célibataires (un pas ascendant est noté a on A selon qu'il est ordinaire 
ou célibataire ; idem pour un pas descendant, noté h ou B), qui séparent les 
équerres, la fonction construit de façon répétitive des équerres, décalées l'une 
par rapport à l'autre de +1 ou de +2, selon que la lettre qui sépare deux 
équerres dénote un pas célibataire ascendant ou descendant : voir la fin de 



la preuve du théorème 11, et la preuve du théorème 14. La variable globale 
compteur^ incrémentée à chaque lecture d'une lettre du mot, permet d'en 
repérer la fin. 

typedef struct {int x, y;} Cellule; 
Cellule Animal [TAILLE_MAX] ; 



static void empiler (i) 

int i; {int h = -2, j = i + DECALAGE, jmax, k; 
for (k = j - 1; k <= j+1; k++) 

if (h < fibre [k] ) {h = fibre [k] ; jmax = k;> 
h++; if (jmax == j) h++; 
Animal [compteur] .x = i; 
Animal [compteur] .y = fibre [j] = h; 



FiG. 19: Code pour empiler une cellule 



La fonction empiler (i), sur la figure 1ë, est élémentaire : elle ajoute une 
cellule sur la fibre i de l'animal, en fonction des tailles des fibres voisines. 
Les cellules sont stockées dans le tableau Animal; la constante DECALAGE sert 
à convertir les abscisses des fibres, qui peuvent être négatives, en indices de 
tableaux. Si la fibre voisine la plus haute, de hauteur /i, est la fibre i elle- 
même, le réseau est triangulaire, et l'ordonnée de la nouvelle cellule vaut 
h + 2, au lieu de /i + 1 en général. 



6.2 Illustrations 



La fonction animal, de la section 3.1, nous a servi a construire les illus- 
trations de cette section. L'un des auteurs a construit un logiciel d'affichage 
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pour X-Window, en utilisant les outils Xview, pour OpenLook. 

Chacune des figures 21, et 22, contient sur sa partie droite un arbre, 
appelé arbre de guingois par les auteurs dans |Q] . Le parcours en profondeur 
de cet arbre, de gauche à droite, donne le mot associé à l'animal, qui se 
trouve en titre en haut de la fenêtre, avec les pas célibataires notés A ou B ; 
cet arbre résume la construction de l'animal à partir du mot. 

Comme l'algorithme de construction d'un animal ajoute des équerres de 
gauche à droite, on peut les lire de droite à gauche, à condition, chaque fois 
qu'on en lit une (sous-arbre issu du point le plus bas de la fibre la plus à 
droite), de l'enlever avant de lire la suivante. 

Les figures 23 et 24 présentent des animaux de grande taille (5000 cel- 
lules), qui sont obtenus sur une station Sun-4 en un temps de l'ordre de la 
seconde, affichage compris. La figure ^ utilise la possibilité de scinder une 
fenêtre en deux : la moitié gauche présente la partie inférieure de l'animal, et 
la moitié droite la partie supérieure, comme on peut le voir à la position des 
"ascenseurs" verticaux. Le logiciel permet d'agrandir à volonté une partie 
de l'animal (délimitée par un rectangle) : la figure 25 illustre ces possibilités. 
Le logiciel afhche aussi, pour chaque animal, sa taille, sa largeur (mesurée 
sur l'axe des abscisses), sa hauteur (mesurée sur l'axe des ordonnées), et le 
nombre de tirages nécessaire pour obtenir le mot aléatoire associé (ce dernier 
renseignement n'est utile que pour les animaux de source ponctuelle). L'util- 
isateur peut modifier la taille avant de demander la génération d'un nouvel 
animal. 

Les animaux montrés ici ont évidemment étés produits de façon aléatoire, 
mais les auteurs ont sélectionné des exemples qui leur ont semblé instructifs 
ou esthétiques ; l'aspect moyen d'un animal est souvent plus simple. 



7 Perspectives 

Cet article a présenté un calcul sur les animaux dirigés, qui permet de 
découvrir une structure cachée de ces derniers, symbolisée par les arbres des 



figures |20|, |21j, et |2^. Parmi les problèmes ouverts concernant les animaux 
plans, le principal est l'évaluation de leur hauteur moyenne (appelée aussi 
longueur). Plus généralement, on devrait pouvoir étudier plus finement la 
structure de ces animaux, et obtenir par exemple leur dimension fractale, 
ou le nombre moyen de composantes connexes du complémentaire. 

Le problème majeur reste celui des modèles de gaz sur des réseaux en 
dimension d > 1, qui correspondent à certains animaux dirigés en dimension 
d + 1 ; le cas du réseau hexagonal a été résolu par R. Baxter (solution du 
modèle des heagones durs, voir Q), et la série génératrice des pyramides 
sur un tel réseau est algébrique, mais une interprétation combinatoire des 
équations obtenues reste à trouver. Le cas du réseau carré (modèle des carrés 
durs) reste entièrement ouvert. 
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FiG. 20: Exemple de construction d'un animal sur réseau carré 
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FiG. 21: Exemple de construction d'un animal sur réseau triangulaire 
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FiG. 22: Exemple de construction d'un animal de source compacte 
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FiG. 23: Exemple d'animal de grande taille sur réseau carré 
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^[z] BBcaacabaccbabbcbBBaaaacbcacabbccaacacbaaabcbccc... ^ 


(Dessiner^ (f^qmnd'w) ( Nouvel animal^ (Quitter") 
X mini -203 X maxi 397 Taille 5000 Largeur 246 
V mini -1 V maxi 648 Tirages fl339 Hauteur 647 


î / h 


1 1 

A. 

1 1 


.1 \i\ IHI 





FiG. 24: Exemple d'animal de grande taille et de source compacte 
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B 



a 



ni ni L. 
E E * 



E E 

K > 



B 



a 



E E 

K > 



FiG. 25: Agrandissements de régions des figures 23 et 24 
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